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Abstract. Two parameters are developed to analyze the CCD images from ground- 
based and/or space telescopes. The first parameter, deduced from the intensity 
profile of the object sharp, is useful to resolve stars and hot pixels. The second 
parameter shape distinguishes the stars from the background cosmic-ray events 
using geometric characteristics defined by its shapes. The parameters are applied to 
a simulated OMC /INTEGRAL image and a HST image. 
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1. Introduction 

The general problem of the fully automatic star detecting algorithms 
on the CCD images is the successful recognition of real star images 
or star-like sources in general from other types of objects and defects. 
These false objects originate from non-ideal characteristics of the CCD 
instrument (noise, bad column and hot pixels) or from the environment 
(nebulae, galaxy or cosmic ray events known as 'cosmics'). The frequent 
requirement is to find real star images on a diffuse background or in 
a crowded star field. Therefore, the detecting algorithm cannot use 
a simple intensity comparing principle, because the intensity of non- 
stellar objects can be comparable to star intensities. 

A number of different methods for the automatic star detection on 
astronomical images has been used during the last two decades. For 
example, the method based on wavelet transforms applied to satel- 
lite ROSAT image processing was published by Damiani et al. (1997). 
The automatic search and classification were implemented in FOCAS 
routine by Valdes et al. (1995). The PC Vista by Trefers and Rich- 
mod (1989) is a package for photometry with a simple algorithm for star 
detection. The well-known program DAOPHOT by Stetson (1987) con- 
tains two parameters sharp and round and a highly robust algorithm 
for star detection. The determination of these parameters originates 
from specific characteristics of the sources on CCD images, but both 
parameters are not sensitive enough to resolve stars and cosmics on 
the CCD images exposed in focus of a space telescope and/or camera 
such as the OMC experiment on-board INTEGRAL (Hermsen and 
Winkler (1997)). 
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2. Distinguishing hot pixels 



The large part of the information describing the detected object is 
included in its intensity profile, but the profile varies widely with 
weather, telescope focus, etc. changes. The intensity distribution cannot 
be described by an analytical function, because it depends on many 
non-controllable factors. Fortunately, the real measured profiles are (in 
the first approximation) close to the two-dimensional Gaussian profile 
both for the ground-based (Moffat (1969)) and space telescopes. 

From this point of view, we can easy distinguish between the hot 
pixel (a single pixel with very different intensity compared to sur- 
rounding pixels) and the real star. The star image with the Gaussian 
profile is usually spread over several pixels. For illustration, the typical 
1-dimensional profiles of the star and of the hot pixel are shown in 
Figures 1 and 2. The maximum intensity of the object ij M0 is at the 
pixel (zoi Jo)- The center of the star is at position 

X] lij ' * lij ' j 

x c — ^ t y° ~ t ' v-U 

where summations are evaluated over a small area around the central 
pixel (io,jo) — for example over 5x5 pixels. This notation is also used 
below. Approximation of the observed intensity by the two-dimensional 
Gaussian profile is the next step of our algorithm 

G(i,j) = Goe-^-^+a-s/c) 2 ]^ 2 + b (2) 

The parameter Go is the estimated maximum of the star profile, b is 
the local sky background. The parameter h is the 'width' of the typical 
star and its numerical value must be known a priori. The stars are 
axis-symmetrical objects according to this formula. Note, that the Go 
parameter is equivalent to, but not quite identical with the Stetson's 
parameter Hi j (integer values are used for x c , y c , i.e. x c = io,y c = jo). 
If we suppose that the statistical distribution of the observed intensity 
has a normal distribution, we can use the least square method 

S = Y J (G(td)-I l , 3 f - min (3) 

for the estimation of the unknown Go and (further unused) 6. 

Further, we can estimate the central intensity Iq from the maximum 
intensity ij M0 by subtracting the local background intensity Bi j from 

liojo ■ 

Iq = liojo ~ Biojo (4) 
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Figure 1. The intensity profile for typical star image. Intensity in arbitrary units, x 
coordinate in pixels. The Gaussian best fit is plotted with parameters: Go = 50.1, 
h = 0.64, x c — 92.0, y c = 76.8. The background intensity is 100 and the observed 
intensity peak Io = 41 satisfies Go > io- 
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Figure 2. The intensity profile for a hot pixel. Intensity in arbitrary units, x co- 
ordinate in pixels. The Gaussian best fit is plotted with parameters: Go = 991.3, 
h = 0.64, x c — 10.0, j/ c = 90.0. The background intensity is 100 and the observed 
intensity peak Io = 1000 satisfies Go < Io- 
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where 

R _ (m)^(»(>J(>) 

lojo - £ WlJ 

(*,j)^(*oJo) 

(wjj are arbitrary weights) and compare its numerical value in param- 
eter sharp: 

sharp = —5-. (5) 
G 

The Figures 1 and 2 indicate the typical values of this parameter. Its 
value is in agreement with the inequality Go > Iq for the star and 
with Go < Io for the hot pixel. We can hence deduce the range of the 
parameter sharp as < sharp < 1 for the star. 



3. Distinguishing cosmic-ray events 

The cosmic ray events frequency on the CCD images exposed at a 
ground-based observatory is low and, moreover, its profiles are narrower 
than the star profiles. The image processing algorithm can then easily 
resolve the images of the stars and those of the cosmic-ray events. But 
CCD images from satellites contain a higher number of the cosmic-ray 
event tracks and the numerical value of the sharp parameter, computed 
for the stars, is close to the numerical value for cosmics. The construc- 
tion of the confidentiality test for its distinguishing with only the sharp 
parameter is not possible. 

The star isophotes of the ideal images have the shape of concentric 
circles, but the cosmics or the bad columns are represented in general 
by ellipses with the major semi-axis at some angle to the x coordinate. 
The mathematical description of an object in image is derived from 
this characteristics. 

The most effective way is the description of the object by its 2nd 
moments 

2 _ Y^hjii-Xc) 2 2 _ Y,kj{j-VcT , fi x 



and 



, _ HIij(i-x c )(j -y c ) 

n-xy — f=tt v'/ 



where the summations are evaluated over the same small area around 
(io, jo)- Now, we can compute the lengths of the object semi-axes as 
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the 'eigenvalues' problem. Its lengths satisfy the equation 

A fox foxy 
foxy A foy 

Resolving of this equations for A gives the sizes of semi-axes as 
l r. 



(h x + h y ) ± yj {h x - h y ) 2 + Ahly 



The numerical values of its parameters (lengths of semi-axes) are dif- 
ferent for elongated object, but approximately the same for stars. With 
respect to their character we define the parameter 



shape = 2 



(hx ~ hy) 2 + 4h% 



xy 



(8) 



The limits of this parameter are: 



— h x — > h y , h xy — > : shape — > 0, the star case 

— h x <-> Ziy, /i^ 7^ : shape — > \h xy j(h x + h y ), the cosmics case 

— h x ^ h y ,h xy -> : |s/iape| -> - h y )/(h x + h y )\ > 0, the 
bad column case, shape parameter is then equivalent to the round 
parameter of Stetson (1987) 

The parameter shape describes the shape of the object. The circular 
objects are interpreted this way as stars and the very elongated object 
as cosmic-ray events or bad columns. Note, that the shape parameter 
represents direct generalization of the round parameter described in 
Stetson (1987). The round parameter is capable of measuring elonga- 
tion only parallel to rows or columns of the image. It was intended to 
identify bad rows or columns, not cosmics. 



4. The simulated image 

In this section we describe the use of the parameters mentioned in 
previous sections to analyze the simulated image of the star field. The 
image (see Fig. 3) is generated as the simulation of the field of view of 
the OMC (Optical Monitoring Camera) prepared for the ESA satellite 
INTEGRAL. The PSF function of the stars in the image is defined as 
the two-dimensional Gaussian profile with 70% of the light energy in 
the central pixel, according to the PSF of the flight optics. This results 
in h = 0.64 for the stars in the simulated image, while their centers 
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Table I. The parameters of objects in the simulated image (legend: No. 1 — hot pixel, 
No. 2 — bad column, No. 3 — cosmics, No. 4 — star). 



No. 


x c 


Vc 


h 


h x 


hy 


foxy 


G 


round 


sharp 


shape 


1 


10.0 


90.0 


1000 


0.17 


0.17 


0.00 


991.4 


0.02 


1.01 


0.06 


2 


5.0 


5.0 


1000 


0.07 


1.41 


0.00 


865.3 


-1.81 


0.91 


1.80 


3 


62.0 


15.0 


199 


1.30 


0.62 


-0.46 


194.2 


0.71 


0.92 


1.20 


4 


92.0 


76.8 


41 


0.69 


0.66 


0.03 


50.1 


0.04 


0.80 


0.10 



are uniformly spread over the unit interval. There are also some other 
sources in the image, namely the bad column, hot pixels and a few 
cosmics (see Fig. 3 and Table I). The image is critically under-sampled. 

The image left-bottom corner has coordinates (0, 0) and right-top 
one (100, 100) in pixels. The bad column is situated at x c = 5. The 
hot pixels occupy positions (10,90), (10,20), (20,10), (20,90) and the 
centers of the three cosmics appear at (10,10), (40,10), (62,15). All 
other sources represent real stars. The intensity of the sky is 100 with 
the Gaussian noise 1.1 (all intensities in arbitrary units). Object No. 4 
(star) and object No. 1 (hot pixel) were used for construction of the 
profiles shown in Fig. 1 and 2. 

The computed characteristics for the selected sources in the sim- 
ulated image are listed in Table I. The center of the object is at 
coordinates x c ,y c , see formulae (1), the intensity peak at pixel (io, jo) is 
related to Iq according to (4). The moments computed from (6) and (7) 
are h x , h y , h xy and the parameter Go is obtained by minimization of the 
sum (3). The last three columns contain the values of the parameters 
sharp, round and shape. 

The values of the numerical differences of the sharp parameter for 
star and hot pixel are small because h is small, but the upper bound 
at 1 separates stars and hot pixels to obviously different groups. The 
values of the round parameter demonstrate small difference between 
cosmics and stars. On the contrary the values of the shape parameter 
of the stars (or hot pixels) and of the cosmics as well as that of the bad 
columns are clearly incompatible. 



5. The HST image 

A wide variety of different kinds of objects appears in the simulated im- 
age, but the real images can also contain others unpredictable defects. 
We test the applicability of the parameters on an example of a real 
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Figure 4- The HST image with stars and cosmics (legend: No. 1 — star, No. 2 
— cosmics, No. 3 — star, No. 4 — nebulosity (?), No. 5 — double star, No. 6 — 
cosmics, No. 7 — star). 
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Table II. The parameters of objects in the HST image (legend: No. 1 — star, No. 2 
— cosmics, No. 3 — star, No. 4 — nebulosity (?), No. 5 — double star, No. 6 — cosmics, 
No. 7 — star). 



No. 


Xc 


2/c 


Io 


h x 


hy 


hxy 


Go 


round 


sharp 


shape 


1 


25.0 


13.8 


406.00 


0.33 


0.54 


0.00 


774.39 


-0.49 


0.51 


0.49 


2 


9.0 


38.0 


724.00 


0.41 


0.95 


-0.26 


717.75 


-0.79 


0.96 


1.11 


3 


63.1 


33.0 


90.00 


0.60 


0.60 


-0.01 


116.62 


0.00 


0.70 


0.04 


4 


77.1 


59.1 


419.00 


0.97 


0.98 


0.02 


502.82 


-0.01 


0.77 


0.04 


5 


83.8 


75.9 


450.00 


0.63 


0.76 


0.31 


723.16 


-0.19 


0.59 


0.92 


6 


48.3 


30.2 


124.00 


0.74 


1.24 


0.72 


877.73 


-0.51 


0.12 


1.54 


7 


72.3 


11.9 


307.00 


0.53 


0.38 


0.00 


1447.14 


0.32 


0.20 


0.32 



image (Fig. 4). We have extracted and zoomed a small sub-window of 
the image from HST. The image has been corrected on photometric 
corrections (bias, flat field, etc.) and exhibits no hot pixels or bad 
columns. We set the value of the h parameter to h = 0.4. 

The list of derived characteristics for a selected objects in this image 
is given in Table II. The meaning of the columns is identical with 
Table I. The objects numbers correspond with Fig. 4. The intensity 
of the sky is 8.3 with the Gaussian noise 1.5 (all in arbitrary units). 

Identifying by 'a visual method' leads to classify the objects No. 1, 
3 and 7 as stars, No. 2 and 6 as cosmics, No. 4 as nebulosity and 
No. 5 is probably a double star. The values of the sharp parameters 
are in the range valid for stars (and hence indicating the absence of hot 
pixels). 'The visual method' classification of stars and cosmics confirms 
the values of the shape parameter — values for cosmics significantly 
exceeds 1. The round parameter has little use for successful recognition 
between stars and cosmics. This example can be described as a typical 
one, because the stars, galaxies and cosmics represent the most frequent 
objects in the HST images. The double star case (No. 5) complicates 
this test and it can easily result in the failure of the whole procedure 
— the shape parameter can be near the limit set in the algorithm and 
we have to set it carefully. 



6. Conclusions 

The two parameters sharp and shape described in this paper are widely 
applicable to the robust detection of the star-like sources in images 
provided by array detectors. The primary reason for deriving the pa- 
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rameters was the development of a reliable automatic algorithm for the 
recognition of stars on CCD images. We can choose following criterion 
for this purpose it the intensity of the central pixel of the examined 
object is above an arbitrarily defined 'level of significance' (threshold) 
and the object lying far from the edge of the frame. Then 

(a) the object is a star for: < sharp < 1 and < shape < 1. 

(b) the object is a cosmics for: shape > 1. 

(c) the object is a bad column for: shape > 1. 

(d) the object is a hot pixel for: sharp > 1. 

This definition provides a strict mathematical criterion for detection 
and classification of the various kinds of objects on CCD images. 

We might specify a range of the sharp and shape parameters for a 
star centered on a pixel as < sharp < 1, < shape < 1. The range 
of sharp for a star centered on the edge between two pixels decrease 
to < sharp < 1/2 and one for a star centered on the corner of four 
pixels decrease to < sharp < 1/4, but the range for shape increase 
to < shape < 2 in this both cases. 
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